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!• MATHBMMIICS. CODING MP APPLICATIONS 

1.1 Introduction 

During the f o^lr week period covered by this report , 700 coded programs 
were run on the time allocated to the Scientific and Engineering Coiaputation 
(SScEC) Group representing work that has been carried on in 40 problems. 

Section 1.2 contains progress reports as submitted by the programmers 
together with an indication of the machine time (and, in some cases, of the 
programming time) expended on each problem. 

The following four problems were completed during this periods #142, A 
Study of Shock Waves; #144, Self -consistent Molecular Orbitals; #176, Connect- 
or Provision in Automatic Telephone Exchanges; #187, Response of a Fuel-flow 
Controller. Various phases of problem 167 (Transient Effects in Distillation) 
have been completed and written up in three Master's theses for the Chemical 
Engineering Department. 

Detailed initial descriptions are provided below for problems #185 (A 
Scale of Turbulence), #186 (Tracking Response Characteristics of the Human 
Operator), #190 (Zeeman and Stark Effect in Positronium) , and #194 (An Augment- 
ed Plane Wave Method As Applied to Sodium). 

Work was initiated during this period on two simulated computers to be used 
in the special summer session course 6.531 (Digital Computers: Business Appli- 
cations). The first computer (a single address computer called SAC) is a modi- 
fication of the 1953 summer session computer (described previously under prob- 
lem #140) involving simulated magnetic tape and omitting the floating point 
programed arithmetic. Progress on SAC will be described under problem #196. 
The second computer is a 3 plus 1 address device designed to simulate an 
actual medium speed drum computer. It will be referred to as TAC and described 
under problem 197. 

Development work continues on the comprehensive system of service routines 
(CS), Teiro new proposed additions to the system are described below under prob- 
lem 100. The first involves the use of a so-called *»director" tape which will 
effect the reading-in of programs in their proper sequence with the information 
necessary for the proper running of each program. Thus it may be possible to have 
such a director tape control without any human intervention a large percentage 
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of the S&EC scheduled machine periods, ^e second prc^osiil involves a simplified 
form for displaying the contents of WMI storage on the oscilloscope, If this 
simplified form is feasible, it -will reduce both the time required for the dis- 
play and the number of registers needed to store the display program. 

During this report period, Digital Computer Laboratory tours were conducted 
for three different groups of visitors, A stHnmary of the size and affiliation 
of each group is given under problem I3I. 

1.2 Programs and Computer Operation 

100. The Comprehensive System of Service Routines ; Best, 88 hours; 
Cjorabelic, 68.5 hours; Demur jian, 40 hours; Denman, .19 hours; 
Frankovich, 68 hours; Helwig, 90 hours; Porter, 2 hours; Saber, 
22 hours; Siegel, 68 hours; WWI, 1218 minutes 

The comprehensive system of service routines provides for the auto- 
matic selection and use of utility programs during computer operation. The 
library of utility programs is permanently recorded on magnetic tape and includes 
conversion programs, binary input programs, equipment checking programs^and post- 
mortem programs. Utility selection is accomplished by using a utility control 
program which has been permanently recorded on the magnetic drum. 

The comprehensive system has been operating successfully during the past 
month. 

The group 11 control program is being modified so that it can be instruct- 
ed to perfoiTO particulaT* functions by means of words appearing on a punched 
paper tape (called a director tape). The director tape mode is selected by 
putting caO in the right insertion register and not pressing the upper activate 
button. The director tape itself is placed in the mechanical reader. 

In the director tape mode every transfer of computer control to the grov^ 
11 control program causes one word on the director tape to be read in. The 
control program then selects a particular utility program on the basis of this 
word and executes this utility program. The director tap® vocabulary has been 
choseil so the director tapes, when printed out, look as much as possible like 
the presently used S&EC perfoiTTiance requests. In particular, the sjmibols 

e, ri, ric, rs, sa # 

fb 100-0-0, fb 1234bi5^ fc 100-0-0* f c i23%5- ■ 
can be typed with their usual meaniiig,' 

A memorandum describing director tape procedure in more detail will be 
issued^when the system is working* 

Work is continuing on a modification to the post-mortem program which wiH 
display post-mortem results on the scope and on the generalized magnetic tape 
recording program described in M-2870, 

Staff 

A program for displaying the contents of CM on the oscilloscope in a simp- 
lified form is undergoing test. Each word displayed consists solely of two 
horizontal rows of points, with the presence of a dot indicating a binary 1 
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and the absence of a dot signifying a binary 0, The upper of the two rows 
shows the word itself j the lower row displays its complement as a check. 

This program promises to be many times faster and shorter than the present 
Group 11 scope post-mortem, which displays the contents of CM as octal fractions 
With the extension of the GSII automatic post-mortem facilities to include the 
scope, the need for an elaborate Group 11 program will b© mach lessened. 

It is projected that this program, to be used only in case of emergency, 
will be recorded on Group 11 in place of the existing scope HI at that time, 
provided that tests shota^ the display to be adequately legible and unambiguous. 
The additional registers made available on Group 11 will be used for needed 
ext«isions of the input program* 

Siegel 

101 C, Optical Properties of Thin Metal Films s Denraan, 26 hours; Loeb, 1 hoan 
WWI, 25 minutes 

A few short runs were carried out to caaplete the work described in 
previous reports. A description of the methods and results of the ccsaputations 
is being prepared for publication, 

106 G, MIT Seismic Pro.iect s Grine, 15 hoursi Simpson, 30 hours; Walsh, 
10 hours; WWI, 3^3 minutes 

Tiie eigenroot i?©utine mentioned in the biweekly of June 13 has been used 
successfully during the past period to determine linear operators which aife 
optimum in a certain sense. These operators, used in conjunction with two 
other programs just written (for computing frequency characteristics of opisrat- 
ors) and an old program, have permitted the completion of a signal-noise sepa- 
ration experiment begun several months ago. In this experiment, performed on 
controlled data, the separation depended on a relatj^vely weak assumption about 
the noise-signal mixture, but compared very favorably with the results of refined 
filter techniques depending on almost complete knowledge of the solution, During 
the experiment further avenues of approach were opened up and are being studied. 

A slightly different attack on the same data is also underway. In this . 
method the multiple filter problem is being considered by examining .frequency 
coherency between seismic trace pairs. There have been canputational delays 
in carrying out this experiment, partly due to ccmputiitional difficulties and 
partly to experimental results requiring procedural modification. 

108 C. An Interpretive Program s WWI, 22 minutes; G, Block 

The individual tapes used in developing this program have been combined 
into one master tape. The only work done during this period was that of making 
a variety of tests on this combined tape. 

120 D. The Aerothermopressor s Gavril, 16 hours; Porter, 3*5 hours; 
WWI, 259 minutes 

During the past month, considerable effort has been directed toward a 
revision of the numerical analysis used in the WWI program for the solution 
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of the first order simultaneous differential equations of the aepoth^ri^^pressor 
process. 

Recently, during computations of the aerothennopressor performance under 
conditions of prescribed Mach Nximber, curious effects were observed which could 
not be e3q)lained on the basis of the physical nature of the process. It was 
suspected that the spurious behavior was the result of accumulated error in 
the numerical solution, and a study of the affect of increment size was accord- 
ingly carried out on a typical case. This study indicated that truncation 
error was indeed a serious effect and that it was necessaiy to improve the num- 
erical analysis which heretofore had been satisfactory for the constant-area 
process. 

The initial revision consisted of a simple forward integration scheme in 

which y J. -, was calculated from backward differences in y' ♦ Four points were 
n^ 1 *^ n 

used with starting values determined by Buler's Method 

The results of tests of this procedure indicated extremely serious instability 
in which oscillation inherent in any numerical solution of this particular set' 
of differential equations was exaggerated by the method itself. After a detail- 
ed study of the results, it was concluded that the procedure was unsatisfactory 
for the following reasons? 

(1) The starting values obtained by Euler's Method contained sufficient 
error ( oscillating) to invalidate the derivatives used in the difference 
equations. 

(2) Once an oscillation had begun, the forward integration procedure was 

unable to compensate for it since differences in y' .^ were not used. Since 

n^ ± 

the differences tacitly involve a polynomial extrapolation of the derivative, 

the error was actually enhanced by the method, 

A significantly more involved procedure is currently being programmed, and 
it is felt that no simple non-iterative method has much chance of success for 
the type of differential equations being treated. The new method consists 
of a forward and successive integration scheme and is essentially Method VII 
of Milne "s Numerical Solution of Differential Equations . This procedure should 
not suffer from the disadvantages listed above. Tentative starting values will 
be determined by Euler's Method but will be adjusted by forward differeiie 
relations until the values contain no oscillation. Values of 7 + -. will be 

computed frcm difference equations involving backward differences in y*„+ i t 

and the calculation repeated until no change in y ^, -, is observed. It is 

hoped that the increased computation per increment can be partly compensated by 
a larger increment size. 

Further computations of aerotherraopressor performance are currently being 
postponed while the new numerical analysis is being developed. 

126 G. Data Reductions WWI, 182 minutes ; D. Ross 

A test run of the Mistake Diagnosis Routine operatir^ on the Data Reduction 
Program was finally successful. Comparison of the results with hand calculations 
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showed several technical programming errors which have now been corrected. 
Another test run is now scheduled. 

An attempt to modify the programs of Problem -l?! slightly to use them as 
function~generating and plotting routines has been beset by logical difficulties 
primarily due to the length of time since the programs were first written. Once 
reestablished, the routines will constitute a simple Fourier synthesizer, 

131. Special Problems (Staff Training. Demonstrations, etc.) ; WWI, 260 minutes 

During the past period a total of I3I people have attended three tours of 
Whirlwind* On June 229 6^ members of the summer course on Operations Research 
were given Flexowriter and computer demonstrations and an informal short talk. 
Nineteen businessmen attending the course on Control Problems of the Executive 
received a similar tour on June 28, Thirty-two foreign students from the MIT 
Foreign Student Projects representing 28 countries, along with representatives 
frcBi Mass, Memorial Hospital, Lincoln Laboratory, the Aeroelastic Structures 
Laboratory J, and the Korwegian-Nether lands Atomic Energy Commission were tour 
guests of the laboratory on July 8, 

X32C, Revision. Extension and Testing of the Subroutine Library Used in Programs 
for Obtaining Data for the Numerically Controlled Milling Machine; routine, 
numerical and logical operations s Porter „ 1.5 hours; Runyon, 50 hours; 
WH^ 78 minutes 

At present p two types of programs are being run. They are concerned with 
milling machine tape checking and airfoil template tape preparation. Some 
difficulty^ evidenced by the appearance of spurious ^rror indications in the 
results, was encountered in checking long tapes. This seemed to be due t© 
excessive loading on the motor of the reader since the errors disappeared 
when a reel which allowed the tape to turn more freely was used* 

An apparently satisfactory inin has been made on the wing template program* 
The portion which was causing trouble was bypassed by using results of a 
previous run, 

1^1 o Matrix Diagonalization and Inversion Subroutine ? Denman, 3 hours; 
Mahoneyj 128 hours; WWI, I58 minutes 

The modification of the subroutine (developed under problem 170) for caldalatiiag 
the inverse, or the square root ©f the inverse, of a real sjnnme trie matrix 
by diagonalization has been tested and is working. It has been found, however, 
more convenient to store the matrix elements in a different order and this change 
will soon be incorporated. 

Some consideration is being given to the use of continued fractions in 
computing functions and a routine for finding In x is now being tested, 

1^2 D, A Study of Shock Waves ; WWI„ 35^ minutes 

This problem was completed during this period. A terminating report will 
appear in Summary Report #38, o 
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i¥^ c* s^i;f~qpnsa.gi^^n'<K ^o3.^g\al^r or^%%^l§' i™:r u minutes 

The program has been rewritten successfully and is now ready to be used. 
The production inins will be described under Problem 201 on the ammonia molecule, 

1^7 C • Energy Bands in Crystals ; WWI, 2311 minutes ; D* Howarth 

During the first week of this four-week period, all the production work which 
it was originally intended to carry out was cc»npleted satisfactorily. This con- 
cemed the calculation of eigenvalues at points of highest symmetry in mcaacnat-um 
space for copper. 

The results greatly exceeded expectations in that three to four figure 
convergence was obtained for all the states with the inclusion of between 4 and 
8 terms in the secular equation. Using a larger number of terms has produced 
five-figure convergence in some cases. This is considerably greater accuracy 
than has been obtained by any other existing method for the calculation of such 
eigenvalues . 

So encouraging were the results that it was decided to calculate also 
certain points of lower symmetry in momentum space; results of significant 
accuracy can be obtained for points of no special symmetry by solution of a 
14x14 secular equation. The number of points being calculated and the fact that 
three separate potentials are being used have made the calculation extremely 
lengthy, and has been made possible only by large amounts of available computer 
time. The entire problem is expected to be concluded during the next bi-weekly 
period. 

The routines used in the problem are be^g modified to permit of simple 
use for other materials, since this appears to be so powerful a method. 

155 "^^ Synoptic Climatology ; Demur jian, 3 hours; Deiaman, 3J^o'u3:'s; 
Blood, 15 hours; WWI, 395 minutes 

Approximately l/l5th of the results have been obtained. Difficulty was 
encountered during this period in determining a satisfactory scale factoring 
arrangement. Scale factors that worked well with the first month's data (an 
unusually cold month) resulted in arithmetic checks when applied to the second : 
month's data (a near-normal month). This has now been overcome. It has not 
been tested on the third month's data (an unusually warm month). If difficulties 
continue it may be advisable to rewrite the program for CSII. 

X59 D. Water Use in a Hydroeleotric.SyBteih ? ,WHI, 785 miaateaj J. DT C; Little 

During the past four w©ekis a considerable nunibfer of intei*©sting results ; 
hAve been obtained. Calculations have been perf ©imed f ©r « vmriety oflmm t o 
dam hydro models and a variety ®f criteria for opfciiiilzation. Operatin^ecuryes 
have been calculated and used ^or running this system under the' conditions of 
the historical river flows. Some exploration has been made of the effect ®f 
relaxing various operating restrictions. Also the mesh size ,©f the calculated 
tables has been varied somewhat to study its, effect on the accuracy. 
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162 C« Determination of Phase Shifts from Experimental Gross-sections; 
Campbell 40 hours; WWI, 45 minutes 

As part of a preliminary program to obtain plots of the sum of the squares 
of the percentage errors in the cross section as a function of each phase sh3^, 
one tape was successfully run giving the error as a function of Sg* the D- /^ 

wave phase shifty when the S and P wave phase shifts have constant values. 
This was done at seven energies for one value of the S wave phase shifts and 
two values of the P wave phase at each energy. 

166 C. Construction and Testing of a Delta-wing Flutter Model ? Porter, 

6*5 hours? WWIj 75 minutes 1 S. Gravitz 

Testing and modification of the completed program is continuing. Prelim- 
inary investigations were made of two proposed lattice network models — the 
first to be machined from a I/8 inch thick aluminum alley plate, and the other 
to be produced from 5/32 inch plate, but otherwise geometrically similar. 

For the I/8 inch plate, using the actual proposed physical dimensions as 
the basis for the stiffness input to the systemg reasonable values for the 
deflection influence coefficients were obtained. It was found that negligible 
changes in the stiffness input to the problem resulted in negligible changes 
in the deflection output, but that arbitrary input changes of substantial 
magnitude can cause much greater corresponding changes in the resultant output. 

For the 5/32 inch thick lattice, a simultaneous program of computation on 
VW and construction and testing of the model is now underway. An inconsistent^ 
in the geometry input to the system was detected and the program is. being mod- 
ified to correct this. Xn the process of pinpointing this error, mn independent 
computation program was performed by use of desk calculators and at every check 
point the preliminary results agreed to within an average of eight significant 
figures with the corresponding data obtained frcaa the Wi program. This close 
agreement tends to indicate that the application of Grout's Method for matrix 
inversion in this part of the program does not yield blcsnfm lap i?»«l^J«ij i^^ 
matrix ill-conditioning. Such ill- conditioning had appeared in the latter 
part of the program where Grout's Method is used to solve a set of simultaneous 
equations. 

Plans for the future include determining the influence coefficients of the 
5/32 inch lattice on 'WW and c ©Bearing the results with experimentally determinea 
values , 

167 D, Transient Effects In Distillation s Jordan, 240 hours; 0*Donnell, 

240 hours; K^ers, 80 hours; Polk, 80 hours; Smith, 80 hours 

This problem did not use any machine time during the period June 14 
through July 11* This report covers the machine time used during the period 
May 17 through June I3, About 1/2 of the machine time used during that period 
was used for the problem of take-off in batch distillation. The r®sainder was 
divided about equally between the problems of equilibration in batch distillation 
and transients in continuous distillation. 
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Ify"ers completed a Master's thesis on the second problem. Time was not 
available for obtaining sufficient data. However, it is expected that additional 
data will be taken in the future with the program which is operating successfully. 

Po3k and Smith completed their Master's theses on the third problem. This 
included considerable data and efforts to correlate it. The work is promising 
and Jordan is continuing it with a Master's thesis this sxararaer. Thus far he 
has been working on further correlation of the data available. 

Considerable effort has been made to correlate the data on take-off in batch 
distillation. As an aid to this effort, a program is being written for Wf which 
will solve the limiting case of batch distillation with no hold-up. 

169. Utilizing a General Purpose Digital Computer in Switching-Circuit Desigi^ : 
WWI9 40 minutes; E, C. Hey 

The rule Xl+X'^+Xf* Xl+X'^-where X and V are single variables 
(i.e. vajriables directly related to only one relay)^Y and 2* are products of 
single variables (none of J being in 2) has been replaced by the more general 

^^^® /product of all 

or any of the i _^ , ^ «- 

va riables of A / \ya riab les of A / 

where X and X' are single variables^ A is a product of single variables^ 



" /product of all \ /product of all or A 

Xl or any of the single)* X«[ any of the single 1+ A«X(f^\A]) + X'Cfg^}) 
\ variables of A / \ variables of A 7 



The runs for timing the reduction program have revealed a reduced expression 
that contains an incorrect term. This will be corrected. 

Other methods for reducing functions have been investigated but none seem 
worthwhile programming in the immediate future. This is not meant to imply that 
the present method is better than those considered more recently^but that those 
considered more recently also have certain seemingly insurmountable difficulties, 
in fact, one of the newer methods allows easier analysis of the degree of minimiz- 
ation obtained. These methods will be listed in the thesis as an aid for future 
studies , 

172 B. Overlap Integrals of Molecular and Crystal Physics ; Gorbato, I60 hours; 
WI, 352 minutes 

Since the last report, the main integral calculation programs have beea 
partially rewritten so as to be more time-efficient. These many routines have 
been combined into one master tape for the convenience of future users. 

Using integral values obtained from the basic subroutines, a tight-binding 
calculation has been carried out for graphite. These results r *t present, rereal 
two features? 1) the wave functions and potentials need not be elaborately 
approximated, and 2) the non-K type secular equations arising in the tight*- 
binding method are sensitive to errors in the overlap matrix elements. The 
latter effect is requiring a theoretical re-examination and extension of the 
original problem formulation, and this work is currently in progress. 
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In the course of using the matrix diagDnalization subroutine in the library 
(LSR) J, it was found that a considerable increase in speed was possible for 
mediura~to high-order matrices by extensive rewriting of the program, A working 
version of this improved routine has been submitted to the LSR. The approxiinate 
improvement in speed is in the ratio of (100+ 28n+in^) to (150+ 52a+ 5n2)» where 
n is the order of the matrix. 

173 « C qurse 6.537 Digital Computer Application Practice s ¥HI, 35 minutes; W» Frank 

The theodolite program was runj it was properly processed through the 
machine but the results were all zero. It was run through once more, stopping 
after one complete cycle and a scope post-mortem was obtained. The program 
is now being de-bugged. 

175 C, Impurity Levels iB Crystals ; Wl^ 411 minutes i G. F. Koster 

The program described previously has been altered and increased in efficiency 
with the help of D^ J, Howarth and F* J. Corbato, This gives a saving in time 
of almost 405^, The inversions of the matrices (H-E) are proceeding ©n a product- 
ion basis and the computation of these inverses necessary for the study of im- 
purity levels in nickel is almost half through, 

176 B, Connector Provision in Automatic Telephone Exchanges ; WI, 317 minutes; 

B, Marr(SJs 
This problem was completed during the present report period, A final 
report will appear in Summary Report No. 38 ^ The res\ilts have been included 
in an MS thesis presented t© the Electrical Engineering Department. 

180 B* Crosse orrelation of Blast furnace Input~Output Data s Denman> 1 hour; 
WWIp 15 minutes p R. Mills 

Progress on the final production program continues 1 the correlation and 
transform programs are both working, and cmly two short control routines remain 
to be checked. Although proceeding slowly^ this work should be completed in 
the near future. 

X81 C* Perturbed Coulodib Wave Functions sH«Paul 1,8 hours; Wl^ 3 minutes 

It has been found necessaiy to generate the wave-function out to r * 15 
in order to find the phase shift and normalization factor needed for c®Biparis®n 
of neutron-deuteron and proton-deuteron scattering at low energies. 

I83 D. Blast Response of Aircraft s Denraan, 1 hour; Sternlight 100 hoxirs; 
WI» 829 minutes 

The production runs for this problem have begimj and during the past period, 
two hundred parameter-sets have been evaluated for aircraft blast response. A 
cubic function encountered in this problem is being evaluated by the Newton- 
Raphson iteration method, and for certain exceptional cases the iteration oscillated 
rather than converged. Since the normal iteration always converges from above, a 
cheek will be incorporated to stop the calculation on a given parameter and go 
on if the iteration starts to hunt from below. The graphs of the computer results 
reveal a very satisfactory^ smooth blast response. 
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184 D. Scattering Electrons from Hydrogen.. Integral Equation s WI, 60 minutes ; 

M» Newstein 
The program to perform a two dimensional integral has been completed and 
checked. Programs to compute various integrands are being checked. Rather than 
ccMpute the integrand independently at each point, advantage is being taken of 
the fact that the integrand varies slowly from point to point. For example, 
instead of computingTX" at x and x+Ax» one uses 

This procedure greatly reduceiS the number of operations required. 

185 B» A Sca^e of Turbulence s Porter, 6 hours; WWI, ^6 mlnxites ; J, Howcroft 

The correlation of eddy velocities or fluctuations between two points with- 
in a turbulent flow will generally vary inversely with the magnitude of the 
distance between the points, ^is offers a method by which eddy size may be 
defined, the length of the eddy being that separation of the points at which the 
correlation beccsnes insignificant. 

Atmospheric turbulence covers a wide range of motions with associated length 
scales varying from large disturbances containing most of the energy of the motion 
to those smallest eddies which dissipate energy as heat by action against viscosity, 
J. W» Tukey states that the minimum resolution of fluctuations is fixed by the 
sampling interval on a stationary stochastic process, ^n upper limit is set 
by the length of the sample. 

Tukey "6 analytic filter was applied to seven sets of wind velocity data. 
These data were simultaneous observations of wind speed at a height of about 
2.5 meters made by four standard Hastings hot thermocouple anemometers with a 
time constant of 0.25 seconds. The anemometers were aligned either parallel or 
normal to the mean wind flow with a separation of 1, 2, and 4 meters between 
individual units on five of the observations, and 3» 6» and 12 meters on the 
remainder. The sampling periods were from 392 to 510 seconds long. TJj© sampling 
interval was one second. 

The analytic filter operates on the sequence 
X., Xg* ••» ^i„i» ^i»^i+i» •••»%« i» % ^ forming the following differences; 

a» First order differences of the j-th probe, 



d 



(1) . 



"" V''*'04 *** 04 T ' ' X *" Xf £>p • . • , I &/ ^ I 



b. Second order differences for the j«th probe » 
c* Third order differences for the j-th probe, 
where [^/Z^] represents the integral part of H/2^. 
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According to Tukey, the sequences of differences as formed above provide 
series in which maximum spectral power is concentrated near periods of 2P seconds 
per cycle, where p is the order of the difference. The filter is only a rough 
one since it lacks sharp cut-off characteristics. In this problem, only the first 
three differences are considered. 

Other investigators have shown that the sums of squares of the differences 
defined above are proportional to the variance and energy at their associated 
frequencies. 

Seven parameter tapes were prepared for presentation to Wiirlwind using the 
observed wind sp^d reading in tens of feet. A program was then prepared in which 
the differences previously defined were formed intermediately without print-out. 
For each difference, for each parameter tape. Whirlwind was called upon t© provide 
the foILowing 30 numbers as its contribution to the solution of the problem; 



fe^ 



i = l 




&/2Li 




> 

J. = 1 


d. .d.. 


6/22) 


^/ik 


fe3 

i = 1 


^ij - \k 



j «1, 2. 3* ^ 



j,K "* X, ^, 3» ^ 



J, K ""X, c.^ jy *+ 



J » k *" X, /&, 3» ^ 



where d.. is the i-th difference in the sequence of differences for the j-th 
probe, *^Thus 63O morabers were printed by Whirlwind. 

With the results from Whirlwind, certain statistical tests were made oiJ the 
data in order to ascertain the degree of validity that the d. ^ *!•© normally 
distributed with mean of zero and with homogeneous variances ^f or synoptic 
observations. This was accomplished by considering each difference for each of 
the seven sampling periods separately as a subject for the following battery of 
tests. First, the homogeneity of variances was tested by using the test descriS^ed 
by Hoel for this purpose. Then an Analysis of Tariance was performed to test 
the hypothesis that tljie means from each of the four probes were essentially the 
same, and the Student s "t" Test was applied to show that this mean was not 
significantly different from zero. Finally, a Chi-square Goodness of Fit test 
was applied to the d^j to determine if they could be considered as having a normal 
distribtition with mean zero and a variance equal to the arithmetic mean of the 
variances of the four individual probes. 
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For each of the differences for each sampling period, the correlation 
coefficients for the six possible combinations of the probes were computed. 
These coefficients were plotted against distance^and the indicated length at 
■which the coefficient becomes sensibly zero noted, ^hus, an eddy with a 
predetermined frequency will have a length, or scale of turbulence, associated 
with it. 

In addition t© the usual correlation coefficient, the d. . were first 
squared (by Whirlwind) and then a figure analogous to the ^ correlation 
coefficient computad and plotted against probe separation. Also, the quantities 
^ |d»^. - d„ 1 were plotted against probe separation, ^he last two plotted 

figures ..will be compared to the plotted correlation coefficients to determine 
the degree of agreement between them. 

186 C. Tracking^ Response Characteristics of the Human Operator ; "WWI, 11 minutes; 

J. Blkind ' ~"~ ' ' ""* ' 
In many control systems (fire control, miasHe control, aircraft flight 
controls , etc) the human operator is an important component of the system. To 
design such systems forcoptimum performance we must know the characteristics for 
all important components including the human operator. 

An experiment has been completed to determine the variability in human 
operator tracking response characteristics. Fourier transforms of the correlation 
functions of stimulus and r^ag^n^e; obtained from one test with one subject have 
been computed by UHI, The power density and cross power density functions thus 
obtained were used to find a quasi-linear transfer function for the subject. 
These results seem promising enough and spectral density functions for other 
tests will be computed as soon as the required correlation functions are avail- 
able. 

187 C. Response of a Fuel-flow Controller s WWI, 57 minutes; M. Merwin 

A system of first order nonlinear ordinary differential equations describ- 
ing the response of a fuel-flow controller has been solved by M. Hej^win of the 
MIT Dynamic Analysis and Control Laboratory on WWI making use of CS II. The 
solution did not prove to be physically realistic so a second problem was then 
solved using the same equations with different constants. 

The Gill method, a variation of the fourth order Runge-Kutta technique, 
was used. The system was the following. 

f * c,y^ Cg 

3^3* -e^kv - c^x^ - c-x « c^^)^ + c„ 



■4 ^8^3 / >J13^ if x^ <c 
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where 

"k = f (2) 
0j = constant 

^2 is • quadratic function of x^^ 

190 D, Zeenian and Stark Effect in Positroiiiimi ; WI, 66 minutes; H. Kendall 

A program to determine the Stark and Zeeman perturbed energy levels of 
the transient electron-positron atom ( positronium) for the n*2 group of states 
is underway ♦ 

Using the modified diaganalization procedure of A. Meckler the energy eigen- 
values, eigenvectors, and reciprocal xifetimes of the decaying states are found 
approximately from the complex non-Hermitian Hamiltonian. 

A variation of the present routine will soon he ready. In addition to the 
above information this will yield the branching ratio: between optical de-excitat- 
ion and direct annihilation for each energy eigenvalue, 

192 D. FrequencT and Phase Spectrom Analysis of SeismoeyaBis: WIj 323 minutes » 

W, P. Walsh^ 50 hours 
Overlapping amplitude and phase spectra of five earthquakes were computed 
and plotted according to "a density plot routine previously described. AH ©f 
the results shiswed travel times of Various reflected and refracted P and S wmvea 
and the dispersivov characteristics of Lofe and Rayleigh waves for Atlantic, 
Pacific, and continental paths. Minimxam group velocities of these waves were ob- 
served for the first time, and they should be instrumental in the determination 
of crystal structure. 

In the next few weeks we intend to continue similar analyses on other 
seismograms in order to gain as much evidence as possible to corroborate claims 
that have already been made, 

193 ^» Eigenvalue problem for Propagation of E, M. Waves s Porter, 3.5 hours; 

wn., 16 minutes; H, B. Dwight 

A tape has been made for one part of the problem that normally takes about 
3 hours with a calculating machine. About 10 test values have been computed with 
the tape, They were correct to 6 to 8 significant figures, so that it is very 
doubtful that there is ai^y error in the tape. The tape is arranged at present 
to handle 15 of these problems at once, and as a large number of solutions are 
needed, these are now being asked for in a routine manner by Lincoln Laboratory, 

A program for a second part of Problem # 193 hsis been typed and submitted. 
No. 193-73-3 « Its length is about equal to that of the first part. It will^ 
as now planned, ultimately be put in series with the first part. 
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194 B. An Augmented Plane Wave Method As Applied to Spditm s WI> 159 misutes 

Maiqy properties ©f solids can be explained succes^fally (conduction, seiai- 
conduction, specific heat, magnetic susceptibility, etc.) if the energy levels 
©f an electron in the solid can be obtained* For simplicity workers have asstamed 
the electrons to move, as a first approximation, independently of one another* 
This approximation is known as the band ©r one-'Olectrofi approximation. Even «©♦ 
solutions are extremely difficult. In the past three methods have been used to 
obtain solutions, each method characterized by its particular approxiination. They 
ares the Bloch method (method of tight-binding ) . the plane-^wave method, and the 
cellular method (method of Wigner-Seitz and Slater ) , U&r& recently the method 
of ortl^Qgonalized plane waves has been applied to the problem of solids. Last 



year the method of augmented planOf^waves wi|s proposed. This method embodies the 



better features of all thfe previous methods. • Howarth ts-ho has recently programmed 
the method on WWI has foutid the method which is simple mathepiatically but 
extremely tedious, ideal for machine calculation. 

I3!l applying the method t© sodium the potential of an electron in sodium 
crystal is represented in a unit cell by a spherically symmetric potential bounded 
by a sphere. At the center of the cell of" this sphere there is joined continuously a 
constant potential which extends throughout the remainder of th.e cell. The 
approximate solutions are taken to be, within the sphere, solutions of the 
radial Schr ©dinger equation for the potential and a particular energy B, say; 
wi-yaout the sphere they are taken to be a plane wave. The radial solution is 
joined continuously to the plane wave. The energy E is then determined as a 
value which minimizes the expectation of energy of the solution. These wave 
functions form a set characterized by s.m. energy E and a propagation vector ©f its 
plane wave part. A set of these wave functions are then combined according t© 
combinations which satisfy a standard non-lambda type secular problem of quantum 
mechanical perturbation theory. The set— and so the order of the secular matrix^ — 
is increased until a succeeding addition does not alter the lowest eigenvalue 
within a fixed amount. This eigenvalue is taken to represent a solution ©f the 
problem, 

Howarth has already programmed this method (Problem 14?) and has applied it 
to copper crystal. Details of his programming can be fcsund in the various bi- 
weekly and quarterly reports. Besides applying the program as it stands to sodiurat 
a modified augmented plane-wave method will be programmed. Because it is basically 
the same method it will utilize Howarth' s program almost totally. The reason f®r 
the modification is that in future applications of the method higher energy levels 
of more intractable solids will be desired and though Howarth has found for the 
lower levels of copper that the method is extremely successful ,the higher levels 
can be obtained as accurately as the lower levels only at tli© expense ©fan increa'g?.- 
ingly larger secular matrix which results in an increasingly larger amount of 
machine time per solution. 

The following problems used ocsaputer time but did mot report t 

107 C, (a) Autocorrelation and (b) Fourier Transform ^ 119 minutes 

Evaluate Integrals 

122 B«, Coulofab Wave Punotions 15 minutes 
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136 Hatrix Equations 

li^ C . Digital Methods of D©tecti23g Signal from Noise 

163 C, Ferrite Phase Shifters in Rectangular Wave GuicJss 

171 C. -Jjnproved Power Estimates 

179 C, Transient Temperature of a Box-Type Beam 

188 G. Effect of Gravity on Relative Water Production 
in Oil Reservoirs 



^' minutes 
25 minutes 
31 minutes 
22 minutes 
22 minutes 
18 minutes 

31 minutes 



197 Three Address Computer 

1»3 Operating Statistics 

l«3l Comuuter Time 

The following indicates the distribution of WWI time allocated 
the S&EC Group during the four week period covered by this report, - 



Programs 


155 hours. 


8 minutes 


Conversion 


11 hours B 


17 minutes 


Magnetic Drum Test 


1 hour 9 


3^ minutes 


Magnetic Tape Test 


1 hour. 


50 minutes 


Scope Calibratioil 


1 hour 9 


22 minutes 


PETR Test 




Z minute^ 


Demonstrations (#131) 


^ hours, 


20 minutes 


Total Time Used 


175 hours f 


;^ minutes 


Total Tijne Assigned 


183 hours 


io ininutes 


Usable Timep Percentage 


95*95^ 




Number of Programs 


700 
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2» COMPTOER EnGIHEERIH& 

(Ko McVicar) 

Computer reliability increased during the past biweekly period. 
This was 5 no doubt, partially because of additional precautionary measures 
being taken by those installing and changing terminal equipment to make 
sure that interference with WWI operation is minimized* 

A program is under way to simplify and standardize maintenance 
and trouble-shooting routines for the computer and terminal equipment* 
An aim of this undertaking is the elimination of the need for highly- 
specialized personnel in the diagnosis of the majority of failure types. 
Efforts are being concentrated on marginal-checking techniques and pro- 
grams which simplify trouble location, fesrice weekly classes are being 
held to increase the dissemination of system logic and circuitry among 
the system technicians^, It is felt that an educational program for the 
terminal equipment is especially desirable since this field has been in 
the past the exelusiTe province of a few. 

2*1 WWI System Operation 

(Lo L. Holmes 5 A. J, Roberts) 

During the past biweekly period the computer dependability was 
excellent o 

We ha¥e renewed our efforts to standardize the circuitry of 
the arithmetic-element control. The step-counter output panel has been 
eliminated. The step counter now uses the bus-driver circiaits in A row^ 
fhe add -to-step-counter circuitry was modified and resulted in a con- 
siderable improvement in the step-counter panel'' s margins, fiming checks 
have been recorded for the division operation, fhe divide -control panel 
will be extensively modified. Included in the changes will be the ex- 
clusion of approximately 11 tubes of the panel's pulse-distribution sec- 
tion. 

A considerable amount of the installation-day work has been 
dedicated to the revision of the marginal-checking logic of the in-out- 
control systemo The work will require about 6 more installatiom days. 
These changes, in addition to the recent introduction of several new lOG 
marginal-ehecking-program routines ^ will siBiplif y the detection and fUialy- 
sis of troubles with this section of the computer* 

(E. H. Gould) 

The timing register on real-time clock has been made more ver- 
satile. Another source of pulses (which must fee synchroniissed with the 
computer "START BELAY** pulses) may be substituted for the 60 pps by feed- 
ing the pulses into Jl-1 on the remote-video-switch panel (R^BP) in E7 
and throwing OH the toggle switch on this panel labeled "EXTRAORDHrART 
IHPDT.'* fhe first two flip-flops of the counter which cannot be read ©ut 
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my be bypassed to give a tiiae count with a l/60-seeond increment or to 

count all of the pulses fed into Jl-1 of the R'VSP„ This is done by 

throwing ON the toggle switch on the RVSP labeled "BT-EASS 2^ HaECOTINTElR«» 
Each switch has an indicator light in test control which is lit when the 
switch is QH<, In normal: operation the switches are OIF and the lights 
are GOTe 

2«2 Terminal Equipment 

2o21 Magnetic Drums (H« L<, Ziegler) 

Changeover of the auxiliary drum to electronic head switching 
for writing is now complete. The change-group delay will soon be re- 
duced to a few microseconds g the exact value to be determined by tests 
now being conduetedo If anyone has reasons for wanting the present 32- 
mseCa delay retained teii|)orarily5 he should contact KJen McTicar or Bob 
Gould immediately o 

Plans for the drum-chassis test setup are essentially complete, 
and the sheet-metal portion of the job should be nearly finished by 9 July. 



3„ ADMIHISTRATIQH AMD PERSOmSL 

New Staff (Jo Co Proctor) 

Arthur G-oddard is working as a DIG Staff member in Group 6345. 

Staff Terminations 

Sylvio Desjardins 
Christopher Laspina 
Hubert Henegar 
Edwin Kopley 
Basil Remis 
Guy Young 

New Non~Staf f (R« A, Osborne) 

Ruth Braid is a new messenger girl at the Barta Buildings 

Transferred Non^-Staf f 

Jean Friberg has been transferred to the Travel Office of 
Division 1<, 

Nora McNeil has been transferred to Buildings and Power • 

Mary Sexton has been transferred to the Director's Office, 
where she is working in the Library* 



DlBtributioa List for S Se E C^ or Group 63k^ 



Internal Distribiition 



Co Aciams 
Bo Ardem 
00 JBest 
Ho Briscoe 
Jo Ackley 

Ho Carr 

Bo GombeXlc 

Jo Gox 

Mo Bemiarjian 

lo Bezuoan 

Co Mafeoney 

Ho BeXXo 
Bo H^aarphrey 
Bo Fellows 
C o Pleiaing 
Jo Frajakovicli 
So "Gill 



Po Bagley 
Fo Heart 
DTo Loeb » ¥ 
Jo Forgie - B 
Bo lelson 



H« Netssaon 
W* Ho Thomas 

Jo W« Forrester 
Ro Bo Eirerett 



Bo Gildea 
Fo Helwlg 
So Hoy 
J, 'KarpczynskL 



Bo Watsom 



Exteraal Blstribrntioa-Blveekly Only 

Bo Bo Ga-vril, 31*^6^^ 

Fo Jo Cor"bat©^ 6-003 

Mo Ho Hellmaa^ Wl~3l6 lastnamentation 

Eo Ho Jacobsen, 6-410 

Jo Bo Co Little, 6-003 



Ro Miller^ 2^-510 

Prof o Morse (T copies), 6-109 
Mareaa (5 Gopjtesi* Ao BoTaiason, 20E-222 
of biweeiay) Bo Boss, Bldgo 32 

No Po Hobbs, 41-219 



Bo Parker 
Jo Porter 
Ho Parechamiaa 
No Saber 
Fo Stwm 
Ao Siegel 
Jo ThompsGH 
¥o Eccles 



Lo Setoit, 4l«-219 

Bo Stemli#t , 33-31^ 
Ho Lin, SS^'Slo 

So Gravitz, 5l-211 

Mo Mo Gben, ii.l*211 

Jo Fo O'BonneH, 12-127 
Bo Marrows, Graduate House 
Bo Campbell, 20B-133 
Mo Co Newstein, 6-003 

Bto Brown, 4-202 
Fo Mo ¥erzuh, 7-306 



So Ho Bodd « B 

Ao Meckler 6-320 

Go Maneen Bldg* B~333 Lincola 

One eopy to Barta Beception Besk 

(Sn-^elope) 

Ao Po Kromer (10 copies of M«Hotes for 



Bldgo 32 « 
Air Force Cambridge Besearcb Center 
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230 Albany Street 

Mathematics Branch 

Office of laval Research 
Washington 25, BoCo 
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